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Nomenclature 


a 

A 

BO 

c 

\ 

F 

Fr 

g 

go 

G 

Gr 

h. 

X 

J 

k 

k 

e 

L 

Nu 

Pr 

q 

Q 

q" 

qo 


acceleration 

area 

Bond number, eqn (7-2) 
specific heat 

blackbody emissive power, eqn (8-2) 
configuration factor, eqn (8-1) 

Froude number, eqn (7-4) 
acceleration of gravity 
standard acceleration of gravity 
irradiation, eqn (8-2) 

Grashof number 
latent heat of fusion 
radiosity, eqn (8-1) 
thermal conductivity 
effective thermal conductivity 

can height or characteristic length (Section 7) 

Nusselt number 

Prandtl number 

heat transfer rate 

heat transfer 

heat flux 

wall heater flux 
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^^> HI 


net 


in 


out 


R 

Ra 

Ra* 


t 

t 

w 

T 



T 

on 


V 

V 
w 
We 
X 
z 
a 


net heat transfer rate into a node, eqn (6-1). 
rate of heat transfer into a node, eqn (6-1) 
rate of heat transfer out of a node, eqn (6-1) 
radial coordinate 
can radius 

Rayleigh number, eqn (7-5) 
modified Rayleigh number, eqn (7-6) 
time 

can wall thickness 
temperature 

average food temperature 

linearization mean temperature, eqn (8-11) 
bulk fluid temperature 
cold-spot temperature 
initial food temperature 

heater deactivation temperature (cut-off) 
heater activation temperature (cut-on) 
volume 
velocity 
thickness 

Weber number, eqn (7-3) 
coordinate, eqn (8-17) 
axial coordinate 
thermal diffusivity 
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a 

absorptivity 


a.o 

thermal diffusivity of water 

P 

volumetric coefficient of 

expansion 

e 

emissivity, eqn (8-2) 

■ 

0 

contact angle, eqn (7-1) 


Mt 

coefficient of dynamic viscosity 

P 

density 


/ 

P 

reflectivity 


a 

surface tension (Section 7) 


Stefan Boltzmann constant 

(Section 8) 

T 

transmissivity, eqn (8-6) 


$ 

size parameter, eqn (9-1) 


Subscripts 


avg 

average value 


ax 

axial direction 


C 

conduction and convection 

effect 

f 

evaluated for the food 


i 

evaluated at radial node i 


j 

evaluated at axial node j 



liquid 


m 

medium of interest 


max 

maximum value 


min 

minimum value 




N 


evaluated at the outer radial node 


R 

rad 

s 

t 

V 

w 

1,2 


radiative effect 

radial direction 

solid 

total 

vapor 

evaluated at the wall^ 
evaluated at surface 1 or 


Superscript 

n evaluated in the n^^ time 


2, respectively 


increment 



6.0 INCLUSION OF PHASE CHANGE (MELTING) IN THE THERMAL MODEL 


There is usually a substantial energy exchange associated 
with a phase change. For example, water evaporates relatively 
slowly from an open, heated, pan because 897.5 Btu of energy 
(at 14,7 psia) must be supplied to each pound of liquid at 212F 
to produce a pound. of vapor at 212F, Hence, considerable energy 
exchange takes place with no temperature increase (at constant 
pressure) . 

Similarly, when ice is required, liquid water is placed in 
a region (e.g., the freezing compartment of a refrigerator) where 
the temperature is well below the freezing point (32F at 14.7 
psia). Even after the temperature of the liquid drops to 32F 
considerable time elapses while the necessary energy is with- 
drawn (approximately 143 Btu/lbm) from the water to cause a phase 
change from liquid to solid. Conversely, when the frozen water 
(ice) is melted, this same energy (143 Btu/lbm), called the la- 
tent heat of fusion, must be supplied in addition to the energy 
which results in the elevation of temperature, called the sensi- 
ble heat, in a single phase. 

6 . 1 The Thermal Model 

The thermal model employed thus far (Section 2.3; Ref: Part 
I) considers only sensible heating of the food and neglects all 
latent heating (i.e., phase changes). For initial temperatures 
above 32F (at 5 psia) and cut-off temperatures below 162F (at 5 
psia), the model is correct. However, when initially frozen 
food is heated, considerable error can be introduced by neglecting 
the amount of energy required to melt the food. This effect was 
first mentioned in Section 4,4 of Part I, To account for the 
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latent heat of fusion the original thermal model described in 
Sections 2 and 3 of Part I was modified. 

In the finite difference solutions, the food was subdivided 
into thirty— six (6— radial and 6— axial subdivisions) toroidal 
nodes (Fig. 4,2; Ref; Part I), (The centerline nodes are actu- 
ally cylinders.) The product of the volume of a given node, its 
density, and its latent heat of fusion is the amount of energy 
required to transform the node from a 32F solid (ice) to a 32F 
liquid. In the case of an initially frozen food the temperature 
of each node is monitored. At the time that the average tempera- 
ture of a given node reaches 32F a record is begun of the net 
energy transferred into the node (i.e., the difference between 
energy entering and leaving) . This net energy can be expressed 
for node (i, j) as: 


^net ^^i,j ^in. . 

If 3 1,3 


- q 


out . 
1/3 


q . = 2rr 

net . , 


1/3 


^ Ar j "\j) 

+ 2rrr.k ,,-2T. ,+t. . 1. 

i Az Li/jfl i/3 i/3~i— i 


( 6 - 1 ) 


The temperature of the node is maintained at 32F until the net 
energy into the node equals the energy required to melt the nu- 
trient. material which the node represents; 

'^5 

2TTr.fizAr p h.^ = ^ At (6-2) 

where Tj is the time increment at which the node first reaches 
32F, Ts is the time increment at which sufficient energy has 
been supplied to melt the node and h^^ is the latent heat 
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(energy) of fusion. 

6.2 Results 

The analysis described in Section 6.1 has been applied to 
the Shylab thermal, food-heating model. Each of the thirty-six 
nodes was monitored. Once the temperature of a given node reached 
32F it was retained at that level until sufficient heat was sup- 
plied to melt it. The temperature then became a dependent vari- 
able in the system again. 

The thermal properties of the food can change substantially 
as a result of a phase change. In particular, as seen in Table 
6-1, the thermal diffusivity decreased by an order of magnitude 
from the solid to the liquid state. The higher value would cor- 
respond to a faster heating rate. In the model, only one (liquid) 
value of the thermal diffusivity is used. This fact could lead 
to a significant discrepancy between the heating times predicted 
by the model and the actual heating times. However, the discre- 
pancy is not so great as it might appear. First, heating time 
in the frozen state is a relatively small portion of the total 
heating time. Also, the food melts first near the heated sur- 
face. Thus the thermal diffusivity quickly drops to the liquid 
value. This liquid partially insulates the rest of the food. 

Heat is not effectively dissipated from the wall region. It 
heats up and cuts off the heater in much the same manner as it 
would even if the proper thermal diffusivity were used for the 
solid state. Hence, the model is conservative in this respect, 
but not significantly. 

Figure 6.1 illustrates the effect of the inclusion of the 
phase change. Both average and cold-spot temperatures are plotted 
against time for both models for the large can. The broken line 
indicates the heater is cycling. The times required to reach an 



TABLE 6-1 


THERMAL PROPERTIES OF SELECTED SKYLAB FOODS [10] 


Food Substance 

Water 

Content 

{%) 

Density 

(lb /ft®) 
m ' 

Latent Heat 
of Fusion 

(Btu/lb ) 
m 

... 

Specific Heat 

(Btu/lb^-F) 

m 

Thermal 

Conductivity 

(BtuAr-ft-F) 

Thermal 
Dif fus ivity 

(ft® /hr) 


it * 

t 

* * 

t 

* * 

t 

Prebuttered Roll 

* 

23.8 

* 

15. 

34. 

0.70 

0.34 

0.19 

0.45 

0.018 

0.088 

Coffee Cake 

* 

25.0 

* 

15. 

36. 

0.70 

0.34 

0.19 

0.46 

0.018 

0.088 

Filet Mignon 

* 

63.2 

58. 

91. 

0.71 

0.40 

0.26 

0.97 

0.0063 

0.042 

Chili with Meat 

* 

66.9 

58. 

96. 

0.74 

0.40 

0.26 

0.99 

0.0061 

0.043 

Prime Rib of Beef 

* 

12 A 

59. 

104. 

0.78 

0.42 

0.27 

1.04 

0.0059 

0.042 

Lobster Newburg 

* 

75.5 

59. 

109. 

o.ao 

0.43 

0.28 

1.09 - 

0.0059 

0.043 

Stewed Tomatoes 

88.4* 

61. 

129. 

0.91 

0.47 

0,30 

1.25 

0.0054 

0.044 


* 


Experimentally determined for actual Skylab foods 
Unfrozen 


Frozen 


A 





TEMPERATURE (F) 



Fig. 6.1: Effect of Phase Change on Temperature Response 
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average temperature of 130F are approximately 1.35 hours (sensi- 
ble heating only) and 2.35 hours (phase change included) - an 
additional 1.0 hours. 

In Section 4.4, it was estimated that an additional 0.39 
hours would be required if it were assumed that the heater ele- 
ment was uncontrolled (on continuously) for the additional time. 
Since about 2^ times as much time is required it can be seen 
that the heater is actually on only about 40% of the additional 
time . 


iz< 
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7.0 LOW-GRAVITY CONSIDERATIONS 

In a low-gravity environment fluids do not necessarily 
settle to the “bottom*' of their containers. Instead they may 
“float" within the container. The loss of direct contact with 
the walls can significantly effect the heat transfer. There- 
fore in low-gravity situations fluid behavior is an important 
consideration in heat transfer analysis. 

Strictly speaking, the term zero gravity does not mean 
that the gravitational force is zero but rather that all net 
external forces acting on the system are zero. External forces 
other than gravity acting on a space vehicle are solar forces 
(pressure) , centrifugal forces (due to vehicle rotation) , thrust 
forces and aerodynamic drag forces. Table 7-1 provides a sum- 
mary of typical forces acting on a spacecraft. Space vehicles 


TABLE 7-1. Typical Accelerations in a Space Vehicle [11]. 


SOURCE 

CONDITIONS 

a/Uo 

1. Solar Pressure 

Low absorptance 

10"^'’ 

2. Centrifugal (vehicle rotation) 

surface 


(a) to maintain vehicle 

Low earth orbit 

10"'" 

parallel to earth surface 

High earth orbit 

10"® 

(b) limit cycle to maintain 

Angular velocity » 

o 

1 

03 

vehicle oriented towards 

0.05 deg/sec; 


sun or star 



3. Thrust for ullage control or 

Thrust = 10 Ibf 

10'* 

from radioisotope propulsion 



4. Aerodynamic drag (C^ == 2) 

Altitude = 100 N mi. 

10'® 


250 N mi. 

10”® 


. 400 N mi. 

o 
1— 1 


13 < 
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experience small accelerations at all times even if only that 
associated with solar pressure alone. 

Knowledge concerning the behavior of fluids in low-gravity 
environments is based in a large part on experimental evidence. 
Low-gravity environments can be achieved by several test techni- 
ques: (a) drop tower, (b) aircraft trajectory, (c) magnetic 

forces and (d) spacecraft. The majority of low-gravity heat 
t^^^sfer data to date has resulted from drop— tower tests (with 
aircraft data a distant second) . 

7 . 1 Liquid/Gas Orientation 

In studying the behavior of fluids in low— gravity environ- 
ments, a knowledge of the predominent forces affecting the fluid 
must be established. In an Earth environment, the gravitational 
forces dominate surface-tension forces, and liquids settle to 
the bottom of containers and gasses collect above. In low— gra- 
vity environments, surface tension forces may dominate. 

The surface tension of a liquid is dependent not only on the 

liquid itself but also on the surrounding environment (e.g., the 

surface tension of water exposed to air is different from the 

surface tension of water exposed to its own vapor or to some other 

liquid like alcohol). Liquids are characterized as either "wetting" 

or "non-wetting" (as illustrated in Fig. 7.1). The contact angle 

0 is expressed in terms of the surface tensions between the liquid 

and vapor the liquid and solid and the vapor and solid 

a by the relation [ 12 ] 
vs 

^ -1 /^vs ^ls\ 

0 . cos ( ^ . ( 7 ^ 1 ) 

vt 

A completely wetting liquid has 0 = O'" , and a completely non- 
wetting liquid has 0 = 180*". Since surface-tension forces are 

14 < 
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WETTING LIQUID NONWETTING LIQUID 

o'" < 6 < 90 ® 90 ® < e < 180 ® 


Fig, 7.1: Contact Angle for Wetting and Non-Wetting Liquids 

molecular in nature and are independent of gravitational fields, 
the contact angle 9 is not changed in low-gravity environments. 
Drop tests performed by Petrash [13] verify this theory; sketches 
showing the configuration of "wetting" and "non-wetting" liquids 
of various contact angles are presented for a cylindrical con- 
tainer in Fig. 7.2, for situations where the container is 50% 
and 90% filled with liquid in a zero-g environment. Since water 
is a totally "wetting" liquid, it can be reasonably assumed that 
nutrient substances would resemble the characteristics of alcohol 
in Fig. 7.2. It should be pointed out that the gas pocket would 
not necessarily be at the "top", since the concept of "top" and 
"bottom" are meaningless in zero-g. 

A dimensionless parameter relating the gravitational (or 
acceleration) forces to the surface tension forces is the Bond 
number , ^ 



Fig. 7 


2: Liquid/Vapor Interface Configuration in Cylindrical 

Containers During Weightlessness 


I 


ED VAPOR ED LIQUID 

50% FULL 90% FULL 



ALCOHOL (0 = 0") 



TETRABROMOETHANE ( 8 = 40°) 



MERCURY (0=125") 


16< 
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Bo 


pL^a 

a 


(7->2) 


As the Bond number increases above unity, the gravitational (or 
acceleration) forces become more dominant, and the liquid is 
"settled”. As the Bond number decreases below unity, surface- 
tension forces become more dominant. The surface-tension char- 
acteristics then determine the relative position of the liquid 
and vapor in the container. 

Another dimensionless parameter involving surface tension 
is the Weber number 


We 


CT 


(7-3) 


which is the ratio of the inertia (or dynamic pressure) forces to 
the surface-tension forces. Again, if the Weber number is greater 
than unity, the flow or inertia forces become dominant. Figure 
7.3 shows graphically the three flow regimes [12]. A Froude 


We 


1000 

100 

10 

I 




.01 


# 

INERTIA DOMINATED 

✓ 

^ ACCELERATION 

OR GRAVITY 
DOMINATED 


CAPILLARY ^ 
DOMINATED m 



.001 


i 


I 


1 


.001 .01 


I 10 100 1000 


Bo 


Fig. 7.3; Flow Regimes 
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number 

number) 


(equal to the square root of the ratio of Weber to Bond 
of unity, i.e.. 


Fr = 


V Bo ~ 


V 




= 1 


(7-4) 


separates the inertia-dominated region from the gravitational (or 
acceleration) -dominated region. The transition between regimes is 
gradual and the selection of a Bond number of unity to indicate 
transition is only approximate. 

Instability in the liquid/vapor interface has been observed 
at various specific values of the Bond number — the critical Bond 
number. Critical Bond numbers appear in the literature ranging 
from 0 to 15 depending on the geometry and the substances involved. 
The critical Bond number (determined from one-g tests) can be 
extrapulated to lower gravity levels (based on the results of Otto 
[12] for cylindrically shaped containers). For a given interface 
(e.g., air and water) and a given critical Bond number, eqn (7-2) 
defines a value for the critical characteristic dimension, L, for 
the system. (This dimension must have been specified in the evalu- 
ation of the critical Bond number.) For a cylindrical container 
this dimension is the diameter. The critical diameter for a 
water/water -vapor interface in a cylinder (a/p = 73 cm^/sec^) at 
various gravitational levels is shown in Fig. 7.4. 


7 - 2 Low- Gravity Heat Transfer 


In a zero-gravity environment, there are no free-convection 
currents as a result of bouyancy. The dominant mode of heat trans- 
fer should be the conduction mechanism, although some fluid motion 
may occur as a result of temperature gradients. The significant 
dimensionless parameter in low-gravity heat transfer is the Ray- 
leigh number Ra which is the product of the Grashof number Gr 
and Prandtl number Pr 


■i8< 
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gg c (T , - T, ) 

^ W D , ^ 

Ra = GrPr = — (7-5) 

The Grashof number is the ratio of buoyant forces to viscous 
forces. The Prandtl number is the ratio of momentum diffusivity 
to thermal diffusivity. High Gashoff numbers indicate significant 
buoyant (free convection) forces, and high Prandtl numbers indi- 
cate significant molecular activity (hence convection) with res- 
pect to thermal conduction. Hence the value of the Rayleigh num- 
ber can be used to indicate the relative importance of convection 
to conduction. 

For a gas in a spherical container at low pressure, Tyler 
[14] found no significant contribution from the convection mecha- 
nism for Rayleigh numbers less than 650. For a fluid layer con- 
fined between two horizontal plates and heated from below, the 
convection is surpressed when the Rayleigh number (based on the 
thickness of the layer w) is below 1700 [15]; this critical value 
was verified experimentally for water as the fluid. As a numeri- 
cal example [15] with w = 0.1 ft (1.2 in.) and a temperature 
difference of lOOF, convection currents are induced in water if 
9/9o ^ ^.2 (lO”®) . This example does not differ substantially 
from the specific cases of interest. If the upper surface layer 
were not bounded by a solid surface but by a free surface, the 
Rayleigh number must be below approximately 1100 for the conduc- 
tion mechanism to dominate. Eckert investigated vertical layers 
of fluid, which may represent a tall cylinder enclosed by two 
plates of height L separated by a distance w; the top and 
bottom were insulated. The regime of pure conduction ended at 
Ra = 500 L/w where the Rayleigh number uses the characteristic 
length w [15]. 

For a surface dissipating a uniform heat flux per unit area 



q’ f a modified Rayleigh number Ra* involving the Nusselt number 
Nu is employed 

gPc p^L^q" 

Ra* = GrPrNu - , ( 7 _ 6 ) 

For heated vertical surfaces, the modified Rayleigh number char- 
acterizes the various flow regimes according to [11]: 

Ra* < 1 Conduction dominated 
10^ < Ra* < 10^ Viscous Flow 

The modified Rayleigh number is shown in Fig. 7.5 as a function 
of the fraction of Earth's gravity g/g© for water at lOOF for 
the special case where q" =2.0 W/in! and L = 1-1/8 in. Since 
the onset of convection currents is directly related to the Ray- 
leigh number, as the magnitude of the gravitational forces are 
reduced, the convection regime is delayed. 

^ • 3 Laboratory Simulation of Reduced Gravity 

Siegel [15] discusses the laboratory simulation of increased 
or decreased gravity fields by rescaling the size of the experi- 
ment. Since the Rayleigh number is the important pareimeter in- 
cluding the gravitational forces, the effects of changes in the 
gravitational field can be simulated by altering the characteristic 
length L. For the isothermal heating surfaces, the Rayleigh 
number depends on the product gL^ ; for the uniform surface heat 
flux, the modified Rayleigh number depends on the product gL^ . 
Therefore, in simulating an increased or decreased gravity en- 
vironment by rescaling the size of Earth based experiment, the 
dimensions of the experimental test apparatus must be multiplied 
by the factor { 9 / 90 )^^ or (g/go)^''* , respectively, in sirau- 
lating reduced gravity, the laboratory experiment would be scaled 


down. 



• "7.5: Onset of Convection Currei 

Rayleigh Number 



s Based on Modified 



WATER 

T=IOOF - 
q” = 2.0 W/in 
L = M/8 in. 
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8.0 ALTERNATE MODEL 

No real surface is perfectly smooth. Even a polished 
surface has microscopic irregularities. When two surfaces 
^re brought together then, the actual contact area is only a 
^taction of the apparent contact area due to the surface rough- 
ness as depicted in Fig. 8.1. 

Conduction heat transfer depends on a conducting medium 
^nd is proportional to the cross-sectional area (Fourier ' s 
law) . Due to the reduced effective area, conduction heat trans- 
fer between two surfaces in contact is reduced from the level 
expected if no interface were present. This phenomenon is usu- 
ally analyzed in terms of an effective "thermal contact resis- 
tance." The reduced cross-sectional area can be viewed as an 
increase in the resistance to the transfer of heat by conduction. 

If the narrow space between the surfaces is evacuated, 
radiation is the only mode of heat transfer available to trans- 
port heat between parts of the surfaces not in contact. if the 
contains a gas, then all nodes of heat transfer (conduc- 
tion, convection, and radiation) are present. However, the 
thermal conductivities of gases are negligible compared to the 
thermal conductivities of most solids. Only a relatively small 
amount of heat is conducted through the gas. The smallness of 
the gap between the surfaces usually restricts convection cur- 
rents severely. (in true zero-g, of course, no currents exist.) 
Hence, convection heat transfer is usually insignificant. There- 
fore, even with a gas present, radiation is usually the dominant 
mode of heat transfer across the gap. (Depending on the percent 
of the actual contact, the conduction directly between the sur- 
faces may still be the single, most significant transport mecha- 
nism.) 
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8.1 Contact Resistance in Skylab Heater 

There are two possible areas where significant contact 
resistance could be encountered in the Skylab configuration; 
First, the contact between the food can and the cavity wall; 
and second, between the food itself and the can wall. Of the 
two, the latter is potentially more critical. The separation 
of the food from the wall due to zero-g effects has been dis- 
cussed in Section 7. Also, the plastic bag containing the 
food within the can could substantially reduce contact area if 
wrinkling should occur. 

While the simple reduction of heat transfer to the food 
due to poor contact is certainly significant, a potentially 
far more important problem, related to contact resistance, is 
present; in the actual tray, the heater-control thermocouple 
is located in the wall of the cavity. Due to poor thermal 
contact between the food and the wall in the vicinity of the 
sensor, the cavity wall may heat significantly faster than 
the rest of the wall which maintains contact with the food. 
(There is no place for the heat to go.) In this case, the 
sensor responds to a temperature which is not representative 
of the system as a whole. The sensor temperature is much 
higher; the heater is turned off prematurely. Since heat is 
not dissipated from the area of the sensor effectively, the 
heater remains off for an extended time. The food heating pro- 
cess can, therefore, be delayed simply because the heater is 
not activated as intended. 
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parameters are simply not known. For example, the actual con- 
tact area may be a function of time. An extensive analysis is 
therefore not justified. The following is a discussion of the 
analysis required to generate a reasonable model. 

The contact resistance between the can and cavity wall is 
neglected. If the food does draw away from the can, the separa- 
tion is over the "ring area" (2TTrjg(Az)) associated with the basic 
discretization scheme discussed in Section 4. (This assumes 
axial symmetry.) The energy is then transported from the outer 
radial node (this is the toroidal node depicted in Fig. 8.2) 
to the adjacent food node by radiation and conduction through 
the nitrogen gas used to fill the can. Since there is no food 
at the wall, the outer radial node corresponds to the wall plus 
the space between the wall and the food. Axial heat transport 
to and from this node is therefore equivalent to axial conduc- 
tion in the wall. 

The radial radiation/conduction transport and the wall 
axial conduction can be put into a form which can be used in 
conjunction with the original model. This is accomplished by 
determining equivalent thermal conductivities and heat capa- 
cities which account for the phenomenon discussed and by using 
them directly in a non-isotropic "conduction" model. 

8.2.1 Equivalent Conductance due to Radiation : The radiation 

exchange between the can wall and the food is approximated by 
the exchange between two parallel walls of emissivities (equal 
to absorptivities) , and Cg . This model is approximate ; 

however, due to the other uncertainit ies and simplicat ions in 
the model, a more complex radiation model is not justified. 

The net radiation exchange between any two surfaces is 
usually expressed in the form: 
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q — JiAiFi-g — (8-1) 

Ri -g 

where J is the radiosity (total radiation which leaves a 
surface per unit time and per unit area) and Fi _2 3 .nd F 2-1 
are the radiation configuration factors (e.g., Fi -3 is the 
fraction of energy leaving surface 1 which reaches surface 
2) . 

The radiosity is therefore the sum of the energy emitted 
and the energy reflected (with no transmitted energy) , or 

J = € + p'g (8-2) 

where e is the emissivity, the blackbody emissive power 

(OT*^ in radiative equilibrium), p the reflectivity, and the 
G the irradiation (total radiation incident upon a surface 
per unit time and per unit area) . 

The configuration factor is a geometry-dependent variable. 
For exchange between large parallel walls Aj -A 3 , Fi ^2 = 

Fs_i =1 (since all the radiation leaving either surface strikes 
the other surface) , and 

Gi = J 2 ; Gs = Ji (8-3) 


Therefore, in thermal equilibrium 


Jl = El °Tt + 

(1-Gi) Js 

(8-4) 

Js = art + 

(1—02 ) 

(8-5) 

where 



p'=l_a'_T = l 

- a' = 1 - e 

( 8 - 6 ) 

(transmissivity, t, is zero, and 

a' = e in thermal 

equili- 


brium by Kirchhoff's law). From eqns (8-4) and (8-5) 

j, = gi art + (l-ci) 63 otS 

1 - (l-ei) (I-03) 

ss< 


(8-7) 
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J 0 T 2 + (l-£a) e, OTT 

^ 1 - (l-Ci) (l-e^) 


( 8 - 8 ) 


Substitution of eqns (8-7) and (8-8) into eqn (8-1) together 
with the assumptions of exchange between parallel walls 


Furthermore , 


111 -3 ^ ojut -it) 

A (parallel _1 ^ 

walls ei tp ” 


T% - Tt = (T2-T1) (Tp+Ti) ( tI+Ti) 


(8-9) 


so that if Ti ef Tg = T 


To - Ti 


A. 


4T*" (Tp-Ti) 


( 8 - 10 ) 


( 8 - 11 ) 


This approximation is commonly called the "linearizing approxi- 
mation" for radiation exchange. For 

0.9 < T/t <1.1 (8-12) 

the error in eqn (8—11) is less than 12%, For the initially 
frozen food (largest temperature range) T is chosen as 

A _ 150 + 1-10 


T = 


= 80 F = 540 R 


(8-13) 


so that 


450 

540 


.84 < ~ < 1.13 
T 


610 

540 


From eqn (8-9) then 


(8-14) 


or 


4t^ 0 

A (parallel ~ 1 I (Tg-Ti) (8-15) 

walls ^ ^ 


^ (Ts-Ti) 
A (parallel e (Xg-Xj 
walls 


(8-16) 
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where 


k 


40-P (X^-X, ) 

1 1 . 
“ + — - 1 
01 


(8-17) 


is an equivalent conductivity and X represents the coordin- 
ate measured perpendicular to the surface. Equation (8-16) 
casts the radiative heat transfer between parallel walls into 
the form of a "conduction model" (Fourier's law). 


8.2.2 Radial Equivalent Conductance : If true conduction and/or 

convection heat transfer through the medium between the walls 
is significant, these mechanisms can also be put into this same 
form (eqn (8-16) ) . The effects of conduction and convection 
are usually combined into a single effective conductivity (a 
technique for enclosed spaces, e.g. [16]) 



1 -g 
A 


_ V (Tg ) 

[parallel e^ (X 2 -X 1 ) 
walls 


(8-18) 


where k^^ is an equivalent conductivity due to conduction and 
convection and can be found [16] correlated as a function of 
Grashof number for free convection situations. 

In space, with no free convection, only conduction through 
the intervening medium is present and 



(8-19) 


where kj^ is simply the thermal conductivity of the medium. 

The conduction and radiation represent parallel paths for 
heat transfer. The total heat transfer then is 






= k 


(Ta - Ti) 
(Xg - Xi) 


+ k 


(Tg - Ti) 
(Xg - Xi) 
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(Ta 

(Xs 


I±1 

X,) 


( 8 - 20 ) 


^ ^ (Tb - T, ) 
A e (Xs - Xi) 


( 8 - 21 ) 


where 


k = k 
e e 

C 


+ k 


R 


( 8 - 22 ) 


is the equivalent total conductance of the medium. From eqns 
(8-17), (8-19) and (8-22) 


k = 


40T (Xs-Xi) 

1 1 
+ — 1 
ei 


+ k 


m 


(8-23) 


8.2.3 Wall Equivalent Conductance ; The sixteen outer nodes 
(4-side, 5-top, 5-bottom, and 2-corner described in Section 
4.4) represent either food plus wall or nitrogen gas plus 
wall depending upon whether the food has separated from the 
wall (Fig. 8.3). The radial heat transfer from the "side" 
nodes and axial heat transfer into the food from the top and 
bottom nodes is handled with the concept of equivalent con- 
ductivity discussed in Section 8.2.1. The heat transferred 
along the wall is now considered. 

Heat is transferred axially between the wall nodes (N,j) 
and (N,j+1) (Fig. 8.4) by two parallel paths — along the con- 
tainer wall and through the food (or gas if the food has 
separated) . The wall-plus-nitrogen or wall-plus-food routes 
therefore form parallel paths for heat conduction. The heat 
conducted along the two paths between the two side nodes (N,j) 
and (N,j+1) is 




= A k 
w w 


T - T 

N, i+1 N,i 


j+1 


- z 


(8-24) 
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Fig. 8.4: Model for Axial Conduction 
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q 


f 


T 


A. 

f f 


N- j + J- 


j + 1 



(8-25) 


where A and A^ are the cross-sectional areas for axial 
w f 

heat transfer in the wall and food (or nitrogen) , respectively, 

and k and k^ are the thermal conductivities of the wall 
w f 

and food (or nitrogen), respectively. From the electrical ana- 
logy for heat conduction (or from an independent derivation) 
the equivalent conductance for the parallel paths is 


k 


e 

ax 


A k 
w w 


+ A k 
f f 


A 


(8-26) 


where A = A + A^. The total heat transfer between nodes (N,j) 
w f 

and (N,j+1) is then 


q 

A |(N,j) to (N, j+1) 


A k + A _ k _ T„ , 
w w f f N,j+1 


j+1 



A 




+ 

"2 




(8-27) 


Similarly for two adjacent nodes on the bottom of the can 


a 

A |i,0 to i+1,0 


^w 




r . ^ 
i+ 1 



(8-28) 


8.2.4 Equivalent Heat Capacity : In the model developed in 

earlier sections, each node had the same heat capacity. (Actu- 
ally, the thermal diffusivity, a = k/pc was specified. The 
changes in the value of the thermal conductivity, k, have al- 
ready been discussed in Section 8.3. The density, p, and the 
specific heat, c, always appear in this analysis as a product. 
This product is called heat capacity.) Since the previous model 
was homogeneous (modeling the food only) this assumption was 
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valid. However, the current model includes regions of food, 
voids, and can wall. The heat capacity must now be specified 
at each node individually. In particular, the equivalent heat 
capacities of the sixteen boundary nodes must be determined. 
This is easily accomplished by using the volume weighted 
average of the heat capacities of the constituents of the node 


(po)e 


p c V 
m=l m m m 

V 

m=l m 


(8-29) 


where the subscript m indicates the different materials com- 
posing the node. 


8.3 A Constant Wall-Temperature Model 

The model used so far has included a simulation of the on/ 
off temperature control heater. This type of boundary condition 
significantly increased the complexity of the model over that 
used in an earlier study [17]. It was felt, however, that its 
inclusion was necessary. Now, with the concept of equivalent 
thermal properties already developed in this section, it will 
be shown that the added complexity of this heat flux boundary 
condition is unnecessary and that a constant-temperature boun- 
dary condition adequately describes the physical system. This 
simplification results only when the aluminum can wall is in- 
cluded within the outer food node. 

In the following discussion, attention is directed to a 
typical outer node (i.e., one that includes the wall). Rather 
than include the details of the analysis of all possible wall 
nodes (e.g., top, center, corner, side, bottom, etc.), only a 
side-wall node is discussed. Therefore, axial heat transfer is 
parallel to the wall, and radial heat transfer is perpendicular 
to the wall. It is with this limitation (side wall node) that 
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the subscripts "rad" (radial) and "ax" (axial) are used in this 
section. The concepts developed here are general and can, of 
course, be extended to other nodes. Only the details would 
differ. 


8.3.1 The case for the Constant Wall -Temperature Model : The 

aluminum can is formed of rolled aluminum sheet (nominal thick- 
ness 0.01 in.). The cross-sectional area available for axial 
conduction along the wall for the larger can is (see Fig. 8.4) 

A = 2ttr t = 2tt 0.01 = 0.118 in^ 

w w 2 


where 

outer 


t is the wall thickness. For axial conduction in the 
w 

food node, the cross-sectional area is 

A^ ^ 2tt (R - — ) Ar = 3.4 in^ 
f z 


The thermal conductivities for the aluminum wall k and for 

w 

the food k_ are 
f 


k a; 130 Btu/ (ft-hr-F) 
w 

a' k 0.3 Btu/ (ft-hr-F) 

f water 


Then from eqn (8-26) 


k 

e 

ax 


.118 

144 


(130) + 


3 .4 
144 


3 . 5/144 


(0.3) 


4.7 Btu/ (hr-ft-F) 


(8-30) 


If the food separates from the wall, k^ becomes the thermal 
conductivity of the replacement gas. The equivalent conductance 
for axial heat transfer is essentially unaffected since as seen 
in eqn (8-30) , a decrease in the value of 
(ft-hr-F) for nitrogen) causes no significant change in the value 
of kg 
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However, the equivalent conductance for radial heat trans- 
fer is significantly affected if the food separates from the wall. 
With no separation, the equivalent conductance for radial heat 
transfer is simply the thermal conductivity of the food, i.e.. 


k 

^rad 



k ^ a: 0.3 Btu/(ft-hr-F) 
water ' 


(8-31) 


If the food does separate from the wall, the equivalent conduc- 
tance is due to the combined effects of conduction and radiation 
across the gap. if the effect of the plastic inner liner is 
neglected, eqn (8-23) can be used. For 
A 

T s: 600 R 

01 = e e; 0.2 
w 

€o = 1.0 

Ar 

X 2 -X 1 = ^ = 3/16 in. 

k^ = k ^0.02 Btu/(ft-hr-F) 

m nitrogen ' 

k = 4 (.1714) X lQ-° (600)^ \«/12 ^ 

®rad 1/0.2 + 1/1.0 - 1 

a; 0.01 4- 0.02 = 0.03 Btu/ ( f t-hr-F) (8-32) 

It is evident by the relative sizes of the contributions to k<c, 

®rad 

that radiative heat transfer in the gap is about one half as 
effective as conduction in the gas. 

From Fourier's law the equivalent axial and radial heat 
transfer are, respectively. 


q 

q 


rad 


ax 


k 

^rad 


k 

e 

ax 


A ^(AT) 
rad 


^ax 


rad 


ax 


k 

^rad 


2ttr(Az) (AT) 

rad 


e 

ax 


2TTR(Ar) (AT) 

ax 


(8-33) 

(8-34) 
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When eqns (8-30) , (8-31) , (8-32) , (8-33) and (8-34) are combined 


^rad 


^ 0.005 (at) 


no sep 


rad 


^rad 


^ 0.0005 (AT) 


sep 


rad 


ax 


^0.1 (AT) 


ax 


(8-35) 

(8-36) 

(8-37) 


Due to the overwhelming preference for conduction along 
the wall over conduction from the wall (by comparison of eqn 
(8-37) with both eqn (8-35) and eqn (8-36)), it should be ex- 
pected that the wall temperature should remain fairly uniform, 
even over the unheated portion (top). Also, because conduc- 
tion from the wall is so severely inhibited, the wall should 
heat up quickly. 

The time required to heat the can wall to the cut-off 
temperature can be estimated if it is assumed that all the 
heat added initially is used in heating the wall only. 

The voltime of the aluminum wall (large, can) 

V , = 21TR L t + 2 ttr 2 t - 0.35 in^. 
w w w 

For an initial temperature of - 10 F the heat required to raise 
the wall temperature to 155 F is 

Q = p V c (AT) ^1.15 Btu 

The side and bottom of the food is exposed to the heater so that 
the heated area is 

A = 2TTR L + ttr^ ^ 25 in? 
w 

The heater supplies 2 W/in? While the heater is activated, the 
heating rate is 

‘^raax " ^ W/in®. (25 in*.) s 170 BtuAr 
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(8-38) 
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The time required to heat the walls from - 10 F to 150 F is 

“w 

t = ^ 25 sec. (8-39) 

This represents the minimum time required since some of the 
heat would go to the food. 

In summary then, it would be expected that within the 
first minute or so of heating the walls should reach the cut- 
off temperature and remain between T and T for the 

off on 

remainder of the heating process. Therefore, a model using a 

constant wall temperature between T and T should ap- 

off on ^ 

proximate the physical system reasonably well. 

8,3.2 Initial Results : Based on the constant wall-temperature 

model verified in the previous section, a computer simulation 
was developed. Since all the heat transfer has been modeled as 
equivalent conduction, the simulation is essentially that for a 
conduction model with nonisotropic thermal properties. A finite 
difference algorithm based on the thirty-six node discretization 
scheme (introduced in Section 4.4, Part 1) was used. 

The food was assumed to have the thermal properties of 
water C^/Q-o = 1) and the wall temperature was taken as the 
average of the two control temperatures (T^ = 149 F) , Figures 
8.5 through 8.10 depict the time responses of the cold spot 
and the average temperatures for initial temperatures of bO F 
and - 10 F . The first three figures are for the large can; 
the last three, for the small can. 

The major uncertainity in the analysis to this point is the 
degree to which (if at all) the food separates from the wall. 
Because of this uncertainity, simulations were run for three 
different contact areas (between the food and the wall) : (a) 

no separation (100% contact) (Fig. 8.5 and 8.8); (b) partial 
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separation (from outer wall) (Figs. 8.6 and 8.9); and (c) com- 
plete separation (Figs. 8.7 and 8.10). 

8.3.3 Corrected Results ; The results represented in Figs. 8.5 
through 8.10 must now be corrected for two reasons. First, 
due to the constant-temperature boundary condition, the heat 
flux at the wall is not controlled in the model. The results 
obtained must therefore be checked to assure that the wall 
flux required in the model is less than (or in the limit equal 
to) that which can be supplied by the heater. Second, in cases 
where the food separates (totally or partially) from the wall, 
some of the nodes which previously represented food, now repre- 
sent nitrogen gas (since there are still thirty-six total nodes) . 
The amount of food heated in the model is therefore less than 
that actually in the can. The simulation then predicts that 
the food heats more quickly than in the actual system simply 
because the model includes less food than the actual system. 

8 -3. 3.1 Correction for Maximum Wall Heat Flux ; From eqn (8-38) 
the maximum rate of heat addition to the large can is 170 Btu/ 
hr (if the heater is activated continuously) , Therefore, the 
food cannot "heat-up" at a rate which exceeds this value. 

For the initially unfrozen food (T^ = 60 F) , all the heat 
added is "sensible" (no phase change involved) . The maximum 
rate at which the average temperature could increase can be 
found by equating the heat absorbed by the food with the maxi- 
mum dissipation from the heater. in a "small" increment of time. 
At, the average temperature increases by AT. These two quanti- 
ties are related through an energy balance 

pV c (AT) - q (At) 

4S< 


(8-40) 
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where q is the mean heating rate over the time increment. At. 
Therefore, the maximum rate at which the average temperature 
could increase can be approximated by 


/dT\ ^ 

Vdt/max \ At /max pVc 


(8-41) 


where Q is achieved by continuous heating, q was 

max ^ max 

evaluated as 170 Btu/hr for the large can in eqn (8-38) . it 
is 91 Btu/hr for the small can. If the properties of water 
are used for the food, the maximum instantaneous rates of in- 
crease in the average temperature with no phase change are. 


(d^iraax ^ (8-42) 

large can 

(d?)max ^ 470”RAr (8-43) 

small can 


The maximum heating rates are expected in the "no separation" 
cases. The maximum slopes determined in eqns (8-42) and (8-43) 
are included in Figs. 8.5 and 8.8. A comparison of the maxi- 
mum slopes to the time response of the average temperature for 
the initially unfrozen food indicates that in the initial heating 
period (t < 0,05 hr), the model is being heated at a rate which 
exceeds the output of the heater. Therefore, the heating time 
will be extended. However, since this excess heating is indi- 
cated for only a short time, the overall heating time should not 
increase by a significant amount. 

For the initially frozen food, the maximum rate at which 
the average temperature can increase is reduced because a large, 
part of the heat added is required to melt the food. This por- 
tion is called the "latent" heat. Equation (8-40) is therefore 
modified to the form 
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pVc (AT) + pV = q(At) 


(8-44) 


As seen in Fig. 8.5 for the large can, the cold spot finally 
melts after about 0.3 hrs of heating. At this time the aver- 
age temperature has increased to about 95 F, The average 
heating rate required over this initial 0.3 hr period is there- 
fore from eqn (8-44) 


= pVc (105) + pV (143) 
^avg 0 . 3 


= 400 Btu/hr 


But the maximiim heating rate has been established (eqn (8-38) ) 
as 170 Btu/hr, At this lower rate, the minimum time required 


t . 
min 


.3(400) 

170 


0.7 hrs 


to melt the food completely and raise the average temperature to 
95 F is 0.4 hrs more than predicted by the model without a heater 
flux restraint. It is apparent then that the heating curve 
should be displaced to at least 1.7 hrs at T = 95 F. Above 95 F, 
the slope is below the maximum slope without phase change, so 
it is not expected that further displacement of the time coor- 
dinate (beyond 0.4 hrs) would be necessary. A similar analysis 
for the small can indicates that heating time is delayed at 
least 0.25 hrs beyond that predicted in the model 

As pointed out, these time delays represent the minimum 
extension possible for the heating process. It is assumed that 
the heaters remain activated continuously in the initial heating 
phase. This is realistic and the results are expected to be 
valid. However, some further delay is possible if the heater 
is deactivated for any appreciable time in this initial period. 


8. 3. 3. 2 Correction for Reduced Volume in Model ; As indicated 
earlier, separation of the food from the walls in the model was 
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achieved by replacing the appropriate "food" nodes with "nitro- 
gen and wall" nodes. The results depicted in Figs. 8.6, 8.7, 

8.9 and 8.10 therefore indicate the temperature response for 
a food sample which is smaller than the actual one. These 
results can be used, however, to estimate the temperature for 
a food sample of the proper size. This is accomplished by 
assuming that the heating time is increased proportionately with 
the increase of food mass or volume. Therefore, at a given 
average temperature, the heating time predicted by the model is 
^'^Itiplied by the ratio of actual food volume to model food 
volume to determine the corrected time to achieve the given 
average temperature. 

A typical evaluation of this area ratio parameter is pre- 
sented in Fig. 8.11 which depicts the geometry for the complete 
separation of food in the large can. The actual food volume is 

TJR^h =12.5 in^ 

The model volume is 

(R - (L - Az) = 8.05 in^ 

The actual volume is therefore 

12.5 , 

Q Q 5 = 1.55 or 155% 
of the model volume. 

The time for the average temperature to reach 110 F (Fig. 
8.7) is predicted to be 1.61 hr. The volume correction would 
indicate that the time should be increased by 55% or to 2.56 hr 
Similar calcualtions yield volume ratios for the other 
cases described in the figures. The volume ratios are presented 
in Table 8-1. 
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Table 8-1 

Ratio of Actual to Model Volume 

Configuration Volume Ratio 

Large Can - complete separation 1.55 

(Fig. 8.7) 

partial separation 1.15 

(Fig. 8.6) 


Small Can - complete separation 1.61 

(Fig. 8.10) 

partial separation 1.23 

(Fig. 8.9) 


8 • 3.3.3 Results ; While the two corrections discussed in the 
previous sections are only approximate, they are necessary if 
full use is to be made of the model. These two corrections 
(maximum heat flux and volume ratio) were applied to the initial 
results' of Figs. 8.5 through 8.10. The corrected results are 
contained in Figs. 8.12 and 8.13 for the large and small cans, 
respectively. In these figures the time for the initially 
frozen food to reach a given average temperature is plotted 

the percent of food— can contact area, (No separation 
to 100% contact area, and complete isolation of the 
food from the wall corresponds to 0% contact area) . The wall 
temperature is fixed at 149 F and the food is given the thermo- 
physical properties of water. The interesting (and reasonable) 
result demonstrated in the figures is that a reduction in con- 
tact area up to about 50% results in only a moderate increase 
in heating time. Heating times are, however, significantly ex- 
tended as complete separation is approached. This general trend 
should be anticipated if one considers the parallel paths (radia- 
tion and conduction through the gas and direct conduction between 
the food and the wall) available for heat transfer. The direct 
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conduction is far more effective (comparing eqns (8-35) and 
(8-36)). As separation begins, most of the heat is still trans- 
ferred effectively by this path. The decrease in effective 
cross-sectional area available for direct wall to food conduc- 
tion is partially off-set by an increase in the temperature dif- 
ference (food does not heat as fast). However, the incremental 
gain in temperature difference begins to decrease as separation 
continues. As complete separation is approached, the severe 
reduction in the cross-sectional area significantly reduces 
direct conduction heat transfer; the heat is transferred (less 
effectively) along the radiation/conduction path. The overall 
effectiveness of the transfer begins to diminish rapidly until 
(with total separation) the less effective path is "carrying” 
the whole heat load. 

Due to the fact that the cans are "filled" with food, it 
IS very unlikely that complete or even near-complete separation 
would occur. if it is assumed that "considerable" separation 
is unlikely, the constant-temperature boundary condition pre- 
dicts a somewhat shorter heating time than the intermittant flux 
boundary condition of the first model. The intermittant flux 
boundary condition model is therefore judged to be the more con- 
servative model, and it is used in the next section to perform 
some optimization studies. 
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9.0 OPTIMIZATION 

The primary objective of this work was to predict the thermal 
characteristics of a particular configuration. However, the model 
provides the tool to pursue parametric studies. These studies 
are performed by maintaining all but one of the problem parameters 
at their * standard" values while allowing the remaining one to 
vary over a range of values. The standard configuration used 
(except as noted) in the following results is: the large container 

(R = 1-7/8 in., L = l-l/8in.), = 155F, = 143F, qo = 2 W/in? , 

a = 0.005633 ft^/hr and an initial temperature of 60F. 

An important question to be answered by this study was to 
determine if the nutrient material lingered in the critical tem- 
perature zone (45-140F) for more than two hours. The results 
are therefore presented in one of two forms: 1) the time response 

of the average temperature, T, for specific values of a single 
P^^^^^^her (e.g,, heater power level) and 2) the time required to 
achieve a given average temperature for a continuous range of 
values for a given parameter. This second form of presentation 
allows one to determine quickly the total time for which the aver- 
age temperature is within the critical range. For example, for a 
given set of parameters the time required to heat to 90F and to 
140F can be determined directly from the graph. The difference 
in these values is the time interval during which the average 
food temperature was within that range. 

9.1 Effect of Thermophysical Properties 

The only thermophysical property influencing the temperature 
response is the thermal diffusivity, a. (Actually, thermal con- 
ductivity does appear in the heat flux boundary condition, e.g., 

•S4< 
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eqn (4-15) (Ref; Part I) . Since it appears only with , the 
effect of its variation can be considered to be a variation in 
qo . The standard value of thermal conductivity is that of water 
at standard conditions, 0.327 Btu/ft-hr-^F, ) The thermal dif- 
fusivity for most nutrient materials is close to that of water. 
However, as previously discussed the reliability of published 
therraophysical data is questionable. Therefore, instead of 
using nutrients as parameters (e.g., response time of beef), the 
dimensionless ratio o/o-o (where qo is the thermal diffusivity 
of water at standard conditions = 0.005633 ft®/hr) is used. 

Figure 5.1a and 5.1b (Ref: Part I) depict the effect of thermal 
diffusivity on heating time. 

As discussed in Section 6, the thermal model developed does 
not have the capability of changing the value of thermal diffu- 
sivity during a run. Therefore, because the thermal diffusivity 
is much higher in the frozen state than in the liquid state and 
since the liquid value is used, the model heats more slowly than 
the actual food. 

9.2 Effect of Heater Power Level 

If food is heated by a uniform heat source which is not con- 
trolled, the heater would remain activated continuously until the 
desired average temperature of the food had been attained. However, 
if the food near the heating surfaces is to be kept from boiling, 
a control mechanism must be included in the heater circuit. The 
temperature/time relationships for the temperature-controlled 
heater are compared to the uncontrolled heater in Fig^ 5.2a and 
5.2b (Ref: Part I) . Both the average and "coldspot" temperatures 
are depicted as a function of logarithmic time. For the large con— 
tainer, the time required to increase the average temperature to 
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140F is increased from about eleven minutes to over an hour 
when the heater is controlled. 

Figure 9.1 indicates the heating time required to achieve 
a given average temperature as the heater power level is changed. 
As seen in the figure, heating times can be significantly in- 
creased by increasing the heater flux in the controlled heater 
only up to about 0.5 W/in^ Above that level, increasing the flux 
level simply inactivates the heater for longer periods of times 
(Fig. 9.2 for the larger container and Fig. 5.3 for the smaller 
container) and decreases the heating time only marginally. Hence, 
from the consideration of power utilization above, the optimal 
power would suggest a heater flux of around 0.5 W/in^ 

^ ^ Effect of Initial Temperature 

The initial temperature of the nutrient material is the most 
influential factor on heating time requirements. Figure 9.3 de- 
picts the time responses for foods initially in the frozen and 
ambient initial states for a range of thermal diffusivity. As 
expected the initially frozen foods heat slowly. For nutrient 
materials with thermal diffusivity below that of water, the model 
predicts heating times over three hours to reach 140F, However, 
almost an hour is required to reach 45F; so that the total time 
in the critical range (45-140F) is not so long as it first might 
appear . 

^•4 Effect of Container Size 

It is desirable during the heating process to have the 
maximum heat transfer surface area per given volume of food. For 
a cylindrical container heated on the sides and bottom^ this size 
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Fig. 9.4: Time for the Average Temperature to Reach Various Levels as 

the Size Parameter Changes for the Standard Configuration 



Fig. 9.5: Time for the Average Temperature to Reach Various Levels as the 

Heater Control Temperatures Change for the Standard Configuration 



average temperature to various levels is presented as a function 

of T /t ratios where T - T = lOF. Of particular in- 

on off off on 

terest is the fact, increasing the level above 150F/160F does not 
decrease the heating time significantly. 

9 . 6 Discussion 

In this report a thermal analysis of the proposed Skylab 
food heating system has been presented. A finite difference 
model was used to carry out parametric studies to determine the 
effect on heating time of (a) thermal diffusivity, (b) heater 
power level, (c) initial temperature (d) container size, and (e) 
control temperatures- A summary of the results is indicated in 
Fig. 9.6. These curves represent the time required for the stan- 
dard configuration model to be heated to an average temperature 
of 140F as one of the model parameters is allowed to vary. The 
intersection of all the lines represents the Skylab configuration 
for heating water from an initial temperature of 60F. As the 
configuration varies the change in heating time can easily be de- 
termined by following the appropriate curve. In particular, it 
is noted that lowering by one-fourth the heater power level does 
not greatly effect the heating time. 
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Fig. 9.6: Time for the Average Temperature to Reach 14CF as 

Each of the Problem Parameters Vary (one at a Time) 
from the Standard Configuration 
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